Imaging of spinor gases 
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We explore the dielectric properties of spinor condensates. Gases with vector (nematic) order, 
such as spin-1 condensates with ferromagnetic (antiferromagnetic) interactions, display optical ac- 
tivity (birefringence). This optical activity/birefringence is respectively seen for light detuned from 
resonance by a frequency which is not too large compared with the fine/hyperfine structure. Careful 
tuning of the light frequency can isolate either of these two effects. These features can be used to 
image spin textures in clouds of spin 1 bosons, and to discern states showing different types of spin 
order. 
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I. INTRODUCTION 

Superfluids display dramatic rotational properties. For 
example, arrays of quantized vortices have been observed 
in rotating gaseous spin-polarized 87 Rb and 23 Na Bose- 
Einstein condensates (BECs) . The physics of rotating 
systems is even richer in the presence of spin degrees of 
freedom: thanks to the interplay of the spin and trans- 
lational degrees of freedom, simple vortices are replaced 
by more complicated spin textures. These textures have 
been studied in many different contexts, from 3 He-A Q 
to quantum Hall systems |3| and, more recently, trapped 
gaseous BECs 0, H, El I3- They are closely related 
to three dimensional skyrmion and hedgehog structures 
which could be engineered in clouds of trapped atoms @ . 

Theoretical studies have analyzed the structure of spin 
textures in rotating spinor BECs for different values of 
the spin, the trap rotation frequency, and the atom-atom 
interactions Experimental studies, beginning with 

pseudo-spin- 1/2 systems |3 and continuing with spin- 

1 and spin-2 condensates pj, have observed individual 
spin textures. These experiments are complemented by 
studies of spin dynamics in non-rotating spin-1 and spin- 

2 clouds 0. 

In most of these experiments [To| , imaging is performed 
by first separating the different spin components using a 
magnetic field gradient (the method pioneered by Stern 
and Gerlach ^3)- The spin components are then sep- 
arately imaged via conventional optical absorption or 
phase contrast techniques. This approach is intrinsically 
destructive and is unable to probe coherences between 
the different spin components. Moreover, images of the 
individual components are very sensitive to global rota- 
tions of the order parameter 

In the present paper, we propose a different imaging 
scheme which is predicted to provide non-destructive, in- 
situ, images of the spin textures and which directly ad- 
dresses the geometrical features of the BEC order param- 
eter (density, spin, and nematicity profiles) rather than 
its components in some specific basis. The method re- 



lies on the dependence of the dielectric properties of a 
spinor atomic gas on the internal state of the atoms. Fo- 
cussing our attention on the case of spin-1 atoms, we 
derive a simple expression for the local dielectric ten- 
sor eij (x) in terms of the local density p(x) , spin density 
S(x), and quadrupolar (nematic) order parameter iVy(x) 
of the atomic gas. Using this expression for ey(x), we 
then study the propagation of polarized light across an 
atomic cloud: we show how information about the spa- 
tial profile of the spin texture can be retrieved by an- 
alyzing the intensity, phase, and polarization profile of 
the transmitted light. Our technique has analogies in 
other areas of condensed matter physics: the polariza- 
tion of the electrons in a solid-state sample [IJ as well 
as the order parameter of nematic liquid crystals |l3j 
can be measured by looking at polarization changes in a 
transmitted, diffracted, or reflected light. As the order 
parameter of spin-1 atomic systems is a three-component 
complex field, it can simultaneously have non-vanishing 
vector (spin) and tensor (nematic) components: we dis- 
cuss possible ways of isolating the contribution of each of 
them. 



II. THE LIGHT-MATTER COUPLING 
HAMILTONIAN 

The electric-dipole Hamiltonian describing the inter- 
action of light with a two-level atom can be written in a 
second-quantized form as: 



Jdx ^ D ma 0],( x ) Ei(x) il> a (x) 

r)ict 



(i) 



The atomic field operators </>, ; (x) and ?/> Q (x) destroy an 
atom at spatial position x in respectively the rj sublevel 
of the excited state or the a sublevel of the ground state. 
Ei(x.) is the i-component of the electric field operator 
at x and D^a = (e : 77 1 64 \g : a) gives the matrix ele- 
ment of the i-component of the electric dipole moment 
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between the sublevels rj and a of respectively the excited 
and ground state. 

For light frequencies w well detuned from the resonance 
frequency lu of the g — > e transition, we can eliminate 
the excited state and write an effective Hamiltonian for 
the ground state atoms coupled to the electromagnetic 
field: 



Hw 



dx £ U aijf} ( X ) 4t (x) &j (x) ^ (x)j (2) 

aij/3 



where the U tensor is defined as: 

Udjp = T7- — —t D* ia D VJl3 , (3) 
h{uj — u>o) 



the indices a and (3 run over the ground state sublevels 
and rj runs over the excited state sublevels. 

In our case of spin 1 atoms, the matter field is a 
three-component vector and can be written in terms 
of its (operator-valued) Cartesian components as ip = 
(ipxT^yjipz) The local one-body density- matrix 

Tftij (x) = ipj(x)ipj(x.) is the tensor product of the two 
distinct spin 1 objects ip and ijp , so it can be represented 
in terms of three irreducible tensor components with an- 
gular momentum 0, 1 and 2. Their explicit expressions 
in Cartesian components are: 



5^m#(x) = <? + (x) • V>(x) 



(4) 



P(x) 

5j(x) = -i£ jk im k i(x) = -i$(x) X ^(x) (5) 
^( x ) = tH»M x ) +™>ji(x)] - 



P(x) 



(6) 



The totally antisymmetric unit tensor is denoted £jki- 
Physically, the scalar p corresponds to the total density, 
the vector S to the spin density, and the second-rank, 
symmetric tensor N to the nematicity (or quadrupole 
moment). An identical decomposition into the inten- 
sity I, the spin density S and the nematicity J\f can be 
performed for the electromagnetic field one-body density 
matrix defined as /y = EjEj. 

As we assume no external field is present to break ro- 
tational invariance, the Hamiltonian is a scalar (that is 
spin 0) object. Scalar objects can be found in the tensor 
product of two quantities only if the two quantities have 
the same spin. This implies that the Hamiltonian (j2J) can 
be decomposed in the form: 



W,,, : : /</x & J(x)p(x) + 6 1 S(x)-S(x) 

+ & 2 ]TA/- iJ (x)A^(x)], (7) 



where the three real coefficients &o,i,2 depend on the in- 
ternal structure of the atom and on the frequency of the 



light. Comparing J7J) with (J2J, one finds 

b = g^C/ Q /3/3a (8) 

a/3 



(9) 



a/3 



&2 — g ^ {U a p a f3 + U aa j3p) — — ^ U aaaa ■ (10) 



a/3 



In the next section, we give explicit expressions for the 
U tensor and the hi coefficients for a few simple cases. 



III. THE DIELECTRIC TENSOR OF A SPIN 1 
GAS 

The form Q of the light-matter Hamiltonian can be 
used for studying the optical potential induced by light 
on atoms, as well as the refractive index observed by light 
while crossing the atomic sample. 

In the first case, a scalar potential comes from the lo- 
cal intensity 7(x) of light, a pseudo-magnetic field comes 
from the electromagnetic spin S(x), and the A/"(x) ten- 
sor couples to the nematicity N(x) of the atoms. A 
first application of this formalism has been given in |17| , 
where effects arising from a kind of spin-orbit coupling 
induced by the pseudo-magnetic field are analyzed for 
atoms propagating in suitably designed optical lattices. 

In the present paper, we concentrate on the second 
class of phenomena, namely on how the optical properties 
of the atomic cloud depend on the spin state of the atoms; 
in particular we shall discuss how the polarization state 
of a light beam after crossing the atomic sample can be 
used to obtain information on the spin state of the atoms. 



A. The general structure of the dielectric tensor 

Under a mean-field approximation in which the quan- 
tum correlations between the matter and the light field 
are neglected, a wave equation for the eigenmodes of the 
electromagnetic field in a homogeneous medium can be 
directly obtained from the Hamiltonian @: 



(Hu (1 - A) - hckP ± ) E = 



(11) 



where 1 is the identity matrix, the projector Pj_ projects 
orthogonally to the wavevector k and the tensor A has 
been defined as: 



Ajfe — — Uijki (ran) , 



(12) 



where eo (= 1 /47r in cgs units) is the permittivity of free 
space. By comparing this wave equation with the Fresnel 
equation for a generic medium of dielectric tensor e: 



-e - k 2 P ± ) E = 0, 



(13) 
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one obtains the following explicit expression for the di- 
electric tensor of the atomic sample at the frequency ui 
in terms of the coupling tensor U in the Cartesian basis 
and the one-body density matrix m: 



(jk — Sjk Uijkl {ma) . 



(14) 



The same symmetry arguments previously used to 
parametrize the Hamiltonian in the form Q lead to the 
following expression of the dielectric tensor (| 1 41) : 

ejk = 8jk + c (p) 5 jk - ic\ Ejki (Si) + c 2 (Njk), (15) 



where Cj = 

corresponds to the isotropic polarizability of the atoms, 
while the vector term proportional to c\ describes their 
optical activity around the axis defined by the atomic 
spin S. Finally, the tensor term proportional to c 2 gives 
birefringence effects, whose principal axis coincide with 
the ones of the nematicity ellipsoid defined as: 



P 



i. 



(16) 



The matrix elements Oyfe, and hence the coefficients 
Cj can be related to the line width T of the atomic tran- 
sition 01- F° r the fundamental Di 2 transitions of a 
typical alkali atom such as 23 Na or 87 Rb, r w 5 • 10 7 s _1 
and \cj\ « 10 _15 cm 3 T/\u> — loq\. A typical in- situ imag- 
ing setup pj| uses detunings in the GHz range (which are 
large enough that absorption effects can be ignored (2i| ) , 
and atomic densities of order 10 cm -3 , resulting in a 
dielectric tensor which differs from unity by an amount 
of magnitude ||e — 1|| w 10~ 2 . As we shall see in sec- 
tion IIII CI quantum interference effects can drastically 
reduce the magnitude of any given Cj in some specific 
frequency domains. 



c (0) 



i.e. (V, = A /3, c^" = —A /2, and c^ = A with 



Ao = — 2\Dq\ 2 / eoh(uj — ojo). This result has a simple in- 
terpretation: a photon polarized along the Cartesian axis 
j interacts only with the component of the matter field 
along the same axis. 

For F e = 1, the Cartesian (j,k,l = {x,y,z}) compo- 
nents of the D tensor have the form Djki — D\ Sjki and 
therefore one has: 



£jk = Sjk + Ai [pS jk - m kj ] 
= Sjk + A\ 



2 i 

-pOjfe + —Sjkl i>i 



N. 



, (18) 



-2bj/eQ. The scalar term proportional to cq i.e. c ( 



(i) 



2Ax/3, c\ 



(i) 



-Ax/2, and & 2 



with Ax = — 2\Di\ 2 /e h(uj — ui ). In physical terms, this 
means that light interacts only with the matter field po- 
larized along direction orthogonal to the electric field po- 
larization. 

Finally, for F e = 2, the dielectric tensor has the form: 



= A 2 



1 r 1 1 



= A 2 



5 5i 1 

9 pSjk ~ i2 EjklSl + 6 Njk 



(19) 



i.e. 



,(2) _ 



= 5A 2 /9, c[ 2) = 5A 2 /12, and c y 2 > = A 2 /6. The 
coefficient A 2 can be obtained from the electric dipolc 
moment D 2 of the optical transition of the F g = 1 — > 
F e = 2 transition as A 2 = — 2\D 2 \ 2 / eoh(u — ojq). More 
specifically, D 2 is defined as the electric dipole matrix el- 
ement between the ttif — \ sublevel of the ground state 
and the rap = 2 sublevel of the excited state. The matrix 
elements of the transitions between the other sublevels 
are related to D 2 by the appropriate Clebsch-Gordan co- 
efficients. 



(2) _ 



B. The dielectric tensor for a single transition 



Effect of the fine and hyperfine structures of 
the atom 



The coefficients c, in the parametrization l|15fl of the 
dielectric tensor depend on the internal structure of the 
atom. Since the ground state has been assumed to have 
spin F g = 1 , absorbing a photon can bring the atoms into 
excited states of angular momentum F e = 0, 1, 2. In the 
present subsection, we shall derive a simple expression of 
the Ci coefficients in terms of the strength of the optical 
transition for the case where only one excited multiplet 
is involved. For notational simplicity we omit the angu- 
lar brackets denoting expectation values of matter field 
operators. 

For F e — 0, the D tensor in Cartesian components 
— { x iHi z }) nas the form Dojk — DgSjk, so the 
dielectric tensor is: 



tjk = Sjk + Ao rrijk 



An 



= S 



■Ao 



P 1 

3 Sjk + -^Sjki Si + N jk 



(17) 



If more than one excited state is involved in the opti- 
cal process, the effective coupling U of the ground state 
atoms to the light is given by the sum of the contribu- 
tions © of each single excited state, each of them being 
weighted by a factor l/(Wj — cj), where huji is the excita- 
tion energy and lu the light frequency. Clearly, the result 
will be dominated by the states which are closest to reso- 
nance. By tuning to frequencies where the contributions 
from different excited states cancel, one can take advan- 
tage of quantum interference effects to make at least one 
of the coefficients 00,1,2 vanish. 

For example, considerable simplification is found if the 
detuning is large compared to either the fine or hyperfine 
splitting. Consider the fundamental nS — > nP transition 
of an alkali atom such as 23 Na or Rb: the S =1/2 elec- 
tronic spin couples to the excited state's L = 1 electronic 
orbital angular momentum to produce two fine compo- 
nents of total electronic angular momentum J = 1/2 or 
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3/2. The ground state, having instead a vanishing orbital 
angular momentum L = 0, contains a single fine com- 
ponent of total electronic angular momentum J = 1/2. 
Both 23 Na and 87 Rb have a nuclear spin of / = 3/2: 
coupling between the nuclear and electronic spins there- 
fore result in the ground state splitting into two hypcrfinc 
components with F = 2 and F = 1 (in this paper we limit 
ourselves to F = 1). The J = 1/2,3/2 fine components 
of the excited state respectively split into two hyperfine 
components of F — 1,2, and four hyperfine components 
of F = 0, 1,2,3. For these alkali atoms, the fine struc- 
ture separation between the D\ (Si/ 2 — * P1/2) and D 2 
(S1/2 —> P3/2) nnes is °f the order of THz. Hyperfine 
structure, resulting from the much weaker coupling be- 
tween the electron and the nucleus, amounts to a few 
GHz for the electronic ground state and to a fraction of 
GHz for the excited state l21l. 



If the detuning A of the light is much larger than the 
hyperfine splitting Ahf, the nucleus, whose direct cou- 
pling to radiation is extremely weak, is not expected to 
play any role in the dynamics. In this regime, the fre- 
quency denominators in the contributions @ to the U 
tensor coming from the different hyperfine components 
are approximately equal. Consequently, U acts as the 
identity matrix in the space of nuclear states. It im- 
mediately follows that the dielectric tensor (|14fl only de- 
pends on the electronic part of the atomic density matrix, 
TO (e) _ Tr„[m] = y~] m (mi \ m \mi) , where the index mi 
runs over the 21+1 possible nuclear spin states. As the 
total electronic angular momentum of the ground state is 
,/ = 1/2, the decomposition analogous to the ones in Q 
and 11511 now gives two components of angular momen- 
tum respectively and 1, but no component of angular 
momentum 2. In physical terms, this means that, once 
the sum over the hyperfine components of the excited 
state is performed, one has C2 = b 2 = and no birefrin- 
gence effect nor any mechanical coupling of light to the 
nematic order parameter N can be present. The coeffi- 
cients ci and Co will, however, generally be non-zero p^ . 
For a large but finite detuning as compared with the hy- 
perfine splitting, the ratio c 2 /ci t o scales as Ahf/ A. 



If the detuning of the light is also large with respect 
to the fine-structure splitting, then one can neglect the 
coupling between light and the electronic spin. The sum 
in © then traces out all of the spin degrees of freedom 
(both electronic and nuclear) , leaving only the electronic 
orbital angular momentum. Consequently, in addition 
to C2 = b 2 = 0, we have c% = b\ = 0, and the atoms be- 
have as a gas of spherically symmetric scatterers, with an 
isotropic dielectric tensor. In this regime, the mechanical 
coupling between light and the atoms does not depend 
on their initial state. This effect is currently exploited 
in optical traps in order to obtain a confining_ potential 
which traps all spin states in the same way 



IV. THE IMAGING METHOD 

The simple dependence of the dielectric tensor on the 
spin and nematic order parameters l|15[) is an useful start- 
ing point for optically imaging the order parameter of a 
spin 1 atomic sample. As a simple example, consider 
a thin, pancake-shaped, Bose-Einstein condensate which 
is rotating around the symmetry axis z. Depending on 
the relative value of the s-wave scattering lengths ao i2 
in respectively the singlet and quintuplet channels, the 
ground state of the system in the rotating frame show 
completely different textures 5J. As two specific exam- 
ples, figs^, and show the nematicity and spin and 
patterns for the antiferromagnetic case, do < a 2 , and the 
ferromagnetic case, a 2 < ao- In the former case, the 
condensate shows defects such as ^-dislocations in the 
nematic order parameter N, while the spin S vanishes 
outside of the cores of these defects. In the latter case, 
the condensate shows a spin pattern, while the nematic- 
ity ellipsoids are pancake shaped, with their minor axis 
aligned with the local spin. 

We imagine that the imaging beam propagates along 
z, i.e. parallel to the rotation axis, and its polariza- 
tion before interacting with the atoms is given by the 
two-component, generally complex, polarization vector 
Pin — (Px,P y )- We assume the cloud is optically thin 
and that the length scale of the spin pattern is much 
larger than the optical wavelength used for the imaging, 
so that diffraction effects during the propagation through 
the cloud can be neglected and the atomic density matrix 
m can be treated as locally uniform. These assumptions 
are valid for sufficiently detuned light (for typical atomic 
densities it suffices that \w — uiq\ > 100 MHz) and suffi- 
ciently weak traps (who ~ 10 Hz). Smaller scale feature 
can be imaged by allowing the cloud to expand ballisti- 
cally before imaging. 

The polarization of the transmitted beam at the trans- 
verse position x_l after propagation through the cloud is 
given by the following expression in terms of a column 
integral along the line of sight: 



Po«*(x_l) = e 1 ' 



z/c 



1 + — J dz (n(xj 



z) 



(20) 

The assumption that the cloud is optically thin implies 
that the magnitude of the term involving the integral 
is much smaller than unity. The matrix giving the lo- 
cal refractive index at the spatial position (xx, z), is de- 
fined as the square root n(x±,z) = VfN of the re- 
duced dielectric tensor e^ xy ^ in the (xy) plane, defined 
as e^ xy > = P± e P± where P± is the projection operator 
orthogonal to the z axis. As one can easily see from (f 1 51) , 
e (xy) depends on the density p, on the S z component of 
the spin and on the N XX yy X yy X components of the ne- 
matic order parameter only. 

If the incident beam is circularly polarized a± , the den- 
sity p and the spin component S z both only give a phase 
shift. On the other hand, the nematicity can mix the two 
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FIG. 1: Nematic order of a 7r disclination in a rotating spinor BEC with weak antiferromagnetic interactions after free expansion 
for time t. The lengthscale is given by a = d^ a \/l + (^h t) 2 , where dho = \] fi/mwho, and Who is the frequency of the harmonic 
trap in which the disclination formed |2jj| . Typically, a ~ 100 /^m. If Who is small enough, then in-situ imaging (where t = 0) 
is possible ( 87 Rb in an wi 10 = lOs^ 1 trap has dho = 9 /im, and is amenable to in-situ imaging). Panel (a): Two dimensional 
projections of the nematicity ellipsoids 116H are shown on a regular grid. Away from the center the ellipsoids are degenerate 
and appear as rods. Panel (b,c): Images of the disclination using circularly polarized a+ probe light. The intensity of the 
generated <r„ component is shown in (b) and a phase-contrast image after mixing with the incident light (with phase (f> — n/2) 
is shown in (c). Panel (d): Image of the disclination using light linearly polarized along x and a crossed polarizer, when ci = 0. 
Panel (e): The z component of the spin S z , as a function of radial position. Panel (f): Image of the disclination using linearly 
polarized light and a crossed polarizer when the detuning is much greater than the hyperfine splitting so that C2 = 0. 



circular polarizations. We assume that the frequency of 
the imaging light is chosen so to have a significant cou- 
pling C2 to the nematicity N (see Appendix). If one il- 
luminates the cloud with pure er + polarized light of am- 
plitude cx~+ , the a- component of the field after crossing 
the cloud has the amplitude: 

a_ = ^ J dz [N xx - N yy + i(N xy + N yx )}. (21) 

In geometrical terms, a <r_ component is generated as 
soon as the (elliptical) cross-section of the nematicity el- 
lipsoid (|16fl along the z — plane is not circular. Af- 
ter blocking the original er_|_ component with a polarizer, 
the intensity of the transmitted light i"_ = |a-| 2 will be 
proportional to the square of the column integral of the 
reduced nematicity parameter 



ll^ll 2 = (N xx - N yy ) 2 + Re[2Ay 2 



along the xy plane. In fig^D we have plotted the sim- 
ulated intensity profile for a ir dislocation in a trapped 
rotating antiferromagnetic (02 > o,q) spinor condensate. 
In the m z = (1, 0, —1) basis, we have taken a condensate 
wavefunction of the form: 



1 


u(x + iy) 



exp 



2d 2 



(23) 



(22) 



with u = 1.22 and the harmonic oscillator length e?h = 
\J h/mujho is of the order of /xm 0. The correspond- 
ing structure of the nematic order parameter is shown 
in figCJi. In this specific case, the magnitude of the ne- 
maticity ||iV|| is symmetric under rotations around the z 
axis. The principal axes of rotate by 7r upon circling 
the disclination center. 

More information on the structure of the nematic or- 



6 



der parameter can be extracted by using phase-contrast 
techniques. Denoting with < 9 < tt the angle between 
the x axis and the minor axis of the nematicity ellip- 
soid, it follows from (|21[1 that the phase of the generated 
<7_ component makes an angle 7r/2 + 29 with respect to 
the incident er+ component. When this <j_ component 
is mixed with the incident u + component (with a mixing 
phase (f>) in a typical phase-contrast scheme , the devi- 
ation of the detected intensity from its background value 
will be proportional to the quadrature Re[e i<6 a_/a+' ) ]. 
As an example of such an images, we have plotted in 
figEt a simulated image for the </> = ir/2 case. The phase 
of the generated er_ at points x and — x differ by ir as a 
result of the different orientation of the nematic order pa- 
rameter. Changing the mixing phase rotates this image 
by <p. 

In experiments on liquid crystals [l3| . one typically 
images using linearly polarized light. A crossed polarizer 
on the other side of the sample selects out the orthogo- 
nal component of the transmitted light. If the nematic 
order parameter is not parallel or perpendicular to the 
incident polarization plane, the polarization plane is ro- 
tated, and light will be transmitted through the crossed 
polarizer. The same approach can be used here, however, 
the optically active regions where S z ^ can also rotate 
the polarization plane, making it more difficult to ana- 
lyze the images. One can isolate the contribution from 
the nematic order by working at a frequency at which 
the coupling of light to the spin vanishes, i.e. c\ = 0. A 
simulated image is shown in fig^i f° r incident light po- 
larized along x: the detected intensity is maximum along 
the y axis where the nematic order parameter makes an 
angle ±7r/4 with x. 

We can isolate the contribution from the atomic spin 
S by working at a detuning which is large compared to 
the hyperfine structure of the excited state. In this case, 
c 2 = and there is no mixing of the a± components. The 
difference between the phase shifts of the a± components 
results in the rotation of the polarization plane of a lin- 
early polarized beam by an angle 9 TO t , which depends on 
the column integral of S z as, 



a wavefunction in the m z — (1, 0, —1) basis given by: 



w Ci 
~2c 



dz S z 



(24) 



The intensity of linearly polarized light passing through a 
crossed polarizer will be proportional to sin 2 (# rot ), which 
in the limit of an optically thin sample is just 9 2 ot . One 
can thus get a direct measurement of the integrated mag- 
nitude of the z component of the atomic spins. In fig^, 
S z is plotted for the 7r-disclination, and in fig^, a simu- 
lated image is shown. 

Since 9 rot corresponds to the phase shift between the 
a± components, one can also measure it via phase con- 
trast techniques. In fig|2t we show a simulated phase 
contrast image of a texture in a ferromagnetic gas. In 
the weakly interacting limit [{|, this coreless vortex has 



v(x + iy) 2 
tp(x,y) = | u(x + iy) | exp 



\x 2 + y 2 \ 



2d 2 



(25) 



with u = —1.03, v = 0.81, and is illustrated in [2^,b. 
Again, dh w /im is the oscillator length. In the external 
region, the spins points downwards, while they reverse 
direction as one approaches the center. Consequently, 
the image in figEt shows a low intensity region at the 
center, and a higher intensity region on a ring of radius 
1.8a (= 1.8 dho for in-situ imaging). 

If the system shows a repeated pattern of vortices or 
7r-dislocations instead of a single one, the imaging may 
be better performed in the far-field plane. The amplitude 
of the Bragg-diffracted light in the direction (kj_,fe z ) is 
proportional to the Fourier transform of the emerging 
field after interaction with the cloud. As an example, we 
have plotted in figOt the Bragg diffraction pattern for 
the lattice of 7r-disclinations shown in figGJJt; the incident 
light is a + polarized and er_ diffracted light is observed. 
The periodicity of the lattice results in an evenly spaced 
series of isolated diffraction peaks whose strengths and 
geometrical arrangement give the detailed geometry of 
the lattice. 



V. CONCLUSIONS 

In summary, we have presented a novel technique for 
both the in-situ and post-expansion imaging of spin tex- 
tures in gaseous spin-1 atomic Bose-Einstein conden- 
sates. The technique is based on the dependence of the 
dielectric tensor on the local value of the density and of 
the spin and nematic order parameters. By considering 
a series of different spin textures, and detection schemes, 
we demonstrated how a range of physical quantities, such 
as the z component of the spin, and the magnitude and 
the orientation of the nematic order parameter can be 
imaged. 

This polarized imaging technique can also be used to 
distinguish between states with and without spin or- 
der: both nematic and spin singlet states have been in 
fact predicted for antiferromagnetically interacting spin- 
1 atoms in an optical lattice |24j . 
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FIG. 2: Panel (a,b): geometrical structure of coreless vortex in a trapped rotating spinor BEC with weak ferromagnetic 
interactions. Lengthscale a is explained in figure 1. In (a), arrows represent the in-plane component of the spin, while (b) 
shows the radial dependance of the z component. Panel (c): image of the coreless vortex by detecting the phase difference 
between the transmitted o± components. 




FIG. 3: Panel (a): Nematic order in array of n disclinations. The closest distance between disclinations whose cores have the 
same spin polarization is D = y/h/mQ,, where Q is the rotation speed ( 23 Na rotating at = 10s _1 has D — 17/xm). As in 
figure 1, this image scales under free expansion. Panel (b): Diffraction pattern seen in <r_ light when o + light is incident on 
this array of disclinations. k x and k y are the transverse wave vectors of the diffracted light. 



APPENDIX A: QUANTITATIVE REMARKS ON 
ABSORPTION AND THE c 2 COUPLING 

Dispersive imaging techniques, such as this one, re- 
quire that absorption is negligible. Minimizing absorp- 
tion is also crucial if one wishes to make in-situ measure- 
ments, as the absorption and subsequent spontaneous 
emission processes heat the sample. The present imaging 
scheme of the nematic order parameter requires detun- 



ings A = l<j — luq not too large than the hyperfine splitting 
in order to have a significant coupling C2 of light to the 
nematic order parameter. 

Since for the D\ line Ajjf /r is of the order of a few 
tens |U, absorption can be avoided while still keeping 
c 2/co,i of order unity. These detunings are of the same 
order as those used in the experiments of an d the 
signals are therefore readily measured. 
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